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Evaluation of the Energy Distribution Function in Physical Adsorption 
as the Problem of Laplace Transform. I. Evaluation of the Energy 
Distribution Function from the Second Gas—Solid Virial Coefficient 


Wyznaczanie różniczkowego rozkładu energii adsorpcji jako problem transformaty 
Laplace’a. I. Wyznaczanie funkcji dystrybucji z drugiego współczynnika wirialnego 
gaz—ciało stałe 


Определение дифференциального распределения энергии адсорбции как проблемы 
трансформации Лапласа. I. Определение дифференциального распределения 
энергии адсорбции из второго вариального коэффициента газ—твердое тело 


INTRODUCTION 


The investigations of heterogeneity of adsorbents are of extremely 
great importance, for both theoretical and practical reasons. Heterogeneity 
is responsible for many interesting and important phenomena in physical 
adsorption [1, 2]. 

The first problem is to describe quantitatively the adsorbent hetero- 
geneity. When assuming that the adsorbent structure is not perturbed by 
the presence of adsorbate, such a fully quantitative description is simply 
the distribution of adsorption energy in the adsorption space. In other 
words, the function e(r) must be known, where e is the energy of adsorp- 
tion, and r is the position vector. However, such a function is rather 
impossible to be evaluated from the experimental data. 

When mutual interactions between adsorbate molecules may be ne- 


glected, it is sufficient to evaluate only the function 2 A at- av 


дЕ 
where У is the volume of adsorption space or the number of adsorption 


sites, depending upon the adsorption model assumed. It results the from 
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fact that in absence of the mutual interactions, correlation effects in the 
adsorbate phase may be neglected. 

According to our earlier discussion the function y(e) seems to be а good 
quantitative description of adsorbent heterogeneity for the majority of 
adsorption systems investigated in practice. This is probably the reason 
for which the theoretical investigations of heterogeneity, made so far, 
concerned the problem of evaluating the function (s). 

A number of methods have been proposed in order to evaluate (t). 
Apart from calorimetric measurements, all these methods may be divided 
into two following groups: 1) exploration of the temperature-dependence 
of adsorption data, 2) exploration of the pressure-dependence of adsorp- 
tion data. The first possibility is connected with the virial description 
formalism. In the virial treatment of physical adsorption, the number N 
of adsorbed molecules is expanded into activity powers [2—5] 


3 
N= 825Z + 83s Z + В,.2 «P anc (1) 


where z is the absolute activity, and Вазиз is the (|-- 1)-th the gas-solid 
virial coefficient. The virial coefficients are integrals, the integrands of 
which include the energy distribution y(e) [2]. The attempts to evaluate 
y(e) from the temperature dependence of gas-solid virial coefficients will 
be discussed in detail in the theoretical section. It will be shown that the 
evaluation of the function may be reduced to the problem of solving La- 
place transform. 


THEORETICAL DISCUSSION 


In the first part of this paper we shall consider the possibility of 
evaluating y(e) from the temperature-dependence of the second gas-solid 
virial coeficient B», As stated above, only such adsorption system is con- 
sidered, the structure of which may be assumed to be quite inert, and 
temperature and pressure independent. In such adsorption system By, is 
defined as follows: 


ae feo EET e а 


Above e(r) is the potential energy field, provided by adsorbent for one 
single molecule, the center of which is at point r. For practical reasons 
the volume V is divided into two parts: 1) volume V, called the volume 
of the surface phase, where the adsorption takes place; 2) volume Vj, 
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called the volume of the free gas phase, where the adsorption forces are 


assumed to be absent, =(т)=0. With the above notice the integral (2) 
may be rewritten in the form 


bzs "Jew [£m Ја (3) 
5 


The most interesting case is the case of attractive forces. Thus we may 
write 


аб s 


The above equation may be rewritten as follows: 


Арад os 
0 


where == — ғ. 


It follows from the above equation that the expression for the second 
gas-solid virial coefficient is simply Laplace transform of function y(e) 
to function Bs (ИКТ). Function Bs, (1/КТ) can be easily found in expe- 
rimental way, for instance by gas chromatography [6]. 

Thus the problem of evaluating the distribution function is reduced 
here to the problem of inversion of Laplace transform, with the transform 
В,, given in a table form. There exists a number of such methods [7—9]. 
The kind of the method used will depend on the course of experimental 
points. In the numerical calculations the experimental results of Sams 
et al. [10] concerning the adsorption of argon on graphite P33 have been 
used. Graphite P33 is usually considered a highly homogenous adsorbent 
[11]. Therefore its distribution function should be very limited with one 
sharp maximum. The Laplace transform should be then of exponential 
type. For such distribution functions two methods of inversion of Laplace 
transform will be proposed in this paper: 1) an approximate method, 
when one has a small number of experimental points of B5; 2) an accurate 
method, when one has at one's disposal a large number of experimental 
values of B, 


1. Approximate method 
Let us consider first the case of a fully uniform adsorption space Ү,, 


having an adsorption energy equal to £. The distribution function y(s) is 
in this case simply the delta function V,5(e’—«,)’). The transform B», will 
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) 
be then the function А exp( бо ). The nearer to delta functions the 


k 


distribution y(e’) the deviations of B», for the exponential behaviour are 
smaller. For some small deviations from the exponential behaviour we 
may write 


ace) 2 (6) 


27. AY 
A 
where the factor —2. — 1 (correction polynome) describes the small 
A4 АКТ 


deviations from the exponential behaviour. Above ey, А; and А, are mean- 
while unknown parameters, which can be found, when fitting best eq. 
(6) with respect to experimental data for В,,. Now the problem of evalua- 
ting y(e’) is quite simple. It is known from the literature [8] that the 
function given by eq. (6) is Laplace transform of the following distribu- 
tion function 


X (E) = | (7) 
А, exp A6 -&) ] for & «e 


oo 
The above function fulfils two physical requirements that " notes dé’ 
0 


: ! 
= constans and бин х (6 ) =0. This function is not quite correct 
co 


for small values of =’. However, the zone У, exhibiting small values of ғ 
play a rather negligible role in the total adsorption effect. Our approxi- 
mate method is similar to the method of Steele [11], who assumed 
the gaussian distribution in the theoretical description of neon and argon 
on black pearl BP and on graphite P33. If one considers two surfaces 
of areas Spp and Sp33, then the ratio of the values of В,, for these solids is 


2 
925. вР__ _ вр л | "ap (8) 


= —————— ex — 
9.5 paa SP33 P2 кг 
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where Овр is the variance of the gaussian distribution for black pearl BP. 
Since the graphite P33 is assumed to be fully uniform, Up;; is taken to 
be equal to zero. Eq. (8) may be rewritten in the form 


8 5 Ugp |. 
Un 25,8P = 6P + 4 per (9) 
B25 p33 Sp33 


Thus the plot In (Ba, вр/В»з 33p) against (КТ)? should yield a straight 


line with the slope equal to Uso . In such a way one gets the 


variance which is here the parameter describing the heterogeneity. 

Stelle has found that the data for this system, when analyzed in this 
manner, did indeed give a straight line. The variance obtained from the 
slope is quite reasonable. It is concluded that the gaussian distribution 
did seem to be a useful representation for heterogeneity. 

In the case of our approximate method, there are three parameters: 
Ai, А», and ғо, which describe heterogeneity. If the deviations of В». from 
experimental behaviour are small, and may be assumed as linear ones, 
Д A 

А 1 А.К уй 
distribution is given by eq. (7). This distribution is like long-normal 
distribution. 


і.е., тау be described by the straight line then the 


2. Accurate method 


Let us rewrite eq. (5) in the form: 


À са 
Bs = E des (i) | в.) nl 


where B'4(1/kT) is a correction function describing the deviations of В,, 
from the exponential behaviour, due to heterogeneity of the surface. Let 
us consider B'4 a function of the argument (—kT), and let us write it 
in the form of a series 


S п+1 
) ) 
8,5(АТ) =) Op ^^ г) (11) 


Hence, 
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wp п+1 
IT MEER. 


n=0 
where о, = Узал. It is known from the theory of Laplace transforms that 
when Laplace transform is a product of functions, then the transformed 
function is a splite of functions. Let, for instance, fi(s):f.;(s) be some 
Laplace transform. Then the transformed function is equal to 


[^ (0) 5 (¢-t)ae аз 


where tt and f(s), f;(s) are the following Laplace transforms 


5) ^ (8) exp (-st at 


D (14) 
f (s) [& (t)exp (-st)at 
In our case this law leads ‘3 the following results 
£ 
ЖЕ) о) (2-27) ос a 
where : 
X, (t) «d (t- &,) (16) 
and 
2, (t) - ) Sto” (17) 
n=0 


Eq. (17) follows from the following law: It is known from the theory of 
Laplace transforms that if the function f(s) can be expanded into the 


series 
а, 
fs) = Sa (18) 


convergent for |s|>R, then the following expansion is valid 
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oo 
a n 
F(t)= Г € (19) 
n=0 
oo 
n+4 
Next, it is easy to check that the series An, ( M, 7) (—kT)"*1 
n=0 
co 
an ) ) J 
is Laplace transform of the series ni £—£4| . From egs. (15—19) 
n=0 
we get 
) 
^ yi ) " 
apa = an 20) 
x (=) - Че) S (E-r far @ 
е п=0 i 
and hence 
"E n ; 
1 
х. (=) ЭЕ (6 ~E) ё) 
n=0 
со 
п+4 
It follows from the experiment that function a, (-k Г) is quickly 


n=0 

convergent, especiallly at higher temperatures. Since (—kT) is rather 
a large value for moderate and higher temperatures, the coefficients а $ 
must very quickly tend towards zero, with the subscript n increasing. 
This makes it possible to truncate the series in eq. (11) after some power 
of (—kT), and to find the coefficients an, using the best-fit procedure. 
In this way we obtain the coefficients in the expansion for the energy 
distribution function, given by eq. (21). Let us additionally notice that 
the series (21) is even more rapidly convergent than series in eq. (11), 
because of the factor 1/n!. It should be emphasized that eqs. (18) and 
(19) are valid only for some class of functions which fulfil the additional 
requirement that. 


| F(t) | <C exp(clël | Q2) 


14 Annales, sectio AA, vol, XXXI/XXXII 
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This additional requirement is always fulfilled, which follows from the 


fact that beu х(=) =0. 


NUMERICAL RESULTS AND DISCUSSION 


As was mentioned above, the experimental results of Sams et. al. 
[10] will be used here in the calculations of the energy distribution func- 
tion (=). They concern adsorption of argon on highly graphitized carbon 
P33 (2700). These experimental data have been applied for two reasons: 
First, they are the most accurate experimental results for B,, reported 
so far in literature. Second, graphite P33 (2700) is considered to be 
a highly homogeneous adsorbent. Therefore, the possibility of quantita- 
tive detecting its heterogeneity by means of our methods, would be 
а good reference for these methods. The experimental data for В»; were 
rather numerous, and, therefore, it was possible to apply them to the 
accurate method in order to calculate у(г’). 

The calculations were performed on the computer ODRA-1204 in the 
UMCS Department of Numerical Calculations. The results are shown in 
Figure 1. The dashed line denotes the function (e’), obtained by trunca- 


eo 
e 


«S 50 
"Е 
540 Fig. 1. The energy distribution 
$ function for adsorption system 
an argon on graphite P33 2700; 1 — 
the function x(e’) obtained by 
j truncation of the series 
320 оо 
F ) An ) ) n 
g 40 n 1 (=- =.) 
— n=0 
after the fifth term, 2 — the 
500 4000 4500 2000 function x(e’), obtained by trunca- 


The energy of adsorption in cal. /mole tion after the sixth term 
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co 
1 
tion of the series a [* г 7) 2d , after the fifth term, whereas the 
n 
п=0 


solid line denotes (ғ), obtained by truncation after the sixth term. The 
differences seem to be rather small. 

Later, it is seen from Fig. 1, that contrary to the assumption made 
by Steele [11], the graphite P33 is not a highly homogeneous adsorbent. 

We would like to emphasize a basic advantage of the accurate method 
proposed here, in comparison with the methods used by Steele [11], 
Ross and Olivier [1], Hoory and Prausnitz [12], Honig 
[13]. In their investigations of adsorbent heterogeneity, the analytical 
form of y(e') has been assumed a priori, whereas in our accurate method 
the form (ғ) results from the calculations performed. 
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STRESZCZENIE 


Zagadnienie ilościowej charakterystyki niejednorodności adsorbentów przedsta- 
wione zostało w oparciu o wirialny formalizm opisu zjawiska adsorpcji. Wykazano, 
ze problem znalezienia funkcji rozkładu energii adsorpeji z matematycznego punktu 
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widzenia jest problemem rozwiązania transformaty Laplace'a. W pierwszej‘ części 
autorzy podają sposób na wyznaczenie funkcji dystrybucji z drugiego wspólczyn- 
nika wirialnego gaz—ciało stałe. 


РЕЗЮМЕ 


В работе доказано, что с математической точки зрения проблема определе- 
ния дифференциального распределения энергии адсорбции является проблемой 
трансформации Лапласа. Указан новый метод определения дифференциального 
распределения энергии адсорбции. 


